Nanotransformation and current fluctuations in exciton condensate junctions 
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We analyze the nonlinear transport properties of a bilayer exciton condensate that is contacted 
by four metallic leads by calculating the full counting statistics of electron transport for arbitrary 
system parameters. Despite its formal similarity to a superconductor the transport properties of 
the exciton condensate turn out to be completely different. We recover the generic features of 
exciton condensates such as counterpropagating currents driven by excitonic Andreev reflections 
and make predictions for nonlinear transconductance between the layers as well as for the current 
(cross)correlations and generalized Johnson-Nyquist relationships. Finally, we explore the possibility 
of connecting another mesoscopic system (in our case a quantum point contact) to the bottom layer of 
the exciton condensate and show how the excitonic Andreev reflections can be used for transforming 
voltage at the nanoscale. 

PACS numbers: 71.35.Cc,71.35.-y,72.70.+m,73.63.-b 



Transport in electronic bilayer systems has recently re- 
ceived increasing attention due to the possibility of ob- 
serving the formation of quantum macroscopic order in 
these systems. Indeed, when an electron layer and a hole 
layer are separated by an insulating barrier that is suf- 
ficiently thick to prevent inter-layer tunneling but suffi- 
ciently thin to induce interlayer Coulomb interaction, an 
excitonic condensate (EC) is predicted to form>ii^ Such 
condensate is a macroscopic quantum coherent state, in 
which electrons in one layer are bound to move coher- 
ently with holes in the other layer. These predictions 
have been confirmed in several experiments performed on 
GaAs quantum wells separated by an AlGaAs barrier, 
both in the quantum Hall regime at total filling factor 
V = 1, see [3| and, more recently, also at zero magnetic 
field^ 

So far, most theoretical studies on transport properties in 
EC were concentrated on the linear response regime;^ 
with a special focus on Coulomb drag configurations i^ 
Other recent works have considered the case of EC con- 
tacted to superconducting clcctrodesji^ whereas current 
fluctuation properties have only been addressed for sys- 
tems where interlayer Coulomb interaction is present but 
is not strong enough to lead to condensatioui^ 
A remarkable advance in the field of EC is expected 
to arise from graphcnc bilayers. Such ECs are pre- 
dicted to exhibit substantially higher critical tempera- 
ture than ordinary semiconductor realizationS )^i^°'^^ due 
to the weaker screening and the higher electron and hole 
densities that can be achieved in graphene. Quite re- 
cently, systems of two graphene layers separated by a 
thin insulating boron nitride film have been realized;^ 
and transport experiments in these systems may become 
a reality in the near future. 

In this paper we derive the full counting statistics (FCS) 
of an EC bilayer, providing its complete low-frequency 
transport characteristics^^ This enables us to investigate 
not only the nonlinear conductance, but also the current 
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FIG. 1. (color online) Sketch of the experimental setup. The 
double layer EC is contacted with four metallic electrodes. 



noise and the higher current cumulants. To this purpose 
we shall adopt the model developed in Q, and evalu- 
ate the cumulant generating function (CGF) of charge 
transfer via the nonequilibrium Green's function tech- 
niqueJ^ Moreover, we shall take a mesoscopic view on 
drag-counterfiow geometries where the top layer is con- 
tacted by leads at different chemical potentials inducing 
a current in the bottom layer that is also part of another 
circuit)^ In our case we study a quantum point contact 
between the two leads of the bottom layer and explore 
the possibility of transforming current on the nanoscale. 
The system, schematically depicted in Fig. [I] consists 
of an electron-hole bilayer, where each layer is contacted 
to two metallic electrodes. While no inter-layer tunnel- 
ing is assumed to occur, the two layers are coupled via 
Coulomb interaction. The Hamiltonian modeling the sys- 
tem reads H = Hn + Ht + Hec- The term Hn ac- 
counts for the four metallic electrodes, characterized by 
electrochemical potentials ^am Fermi distribution func- 
tions Uaa, and an energy-independent density of states 
po- a = L,R refers to the contacts on the left/right side 
of the bilayer, whereas a = T,B labels the top and bot- 
tom layer, respectively. Ht describes the particle tunnel- 
ing between the layers of the EC and the metallic contacts 

Ht= J2 7a.(4*. + *^a.), (1) 

a-=T, B, a=L, R 

where ^lt,b^ 1rt,b sxc the tunneling amplitudes, 



Lt/b, Rt/b the electron field operators for the four 
leads, and \E'-j, ^ the field operators for electrons in the 
EC layers at the position x — Q,l (for L, R), respectively. 
As spin is irrelevant in the effects we are investigating we 
consider a spinless system. Finally Hec describes the EC 
bilayer. For the EC all important features we want to de- 
scribe are captured by a simple one-dimensional model^ 
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where I is the longitudinal distance between the elec- 
trodes, ^ = (^T, ^b)"^ is the two-layer spinor, Ht (Hb) 
describes the electron (hole) single-particle term of the 
top (bottom) layer. The inter-layer Coulomb interac- 
tion is describediii by an exciton order parameter A(a;), 
which is in general a space-dependent quantity. Its bulk 
absolute value Aq at equilibrium represents the exci- 
tonic gap and determines the excitonic correlation length 
^EC = vfI ^Q- We use units such that fc^ = e = fi, = 1 
and Go = 2eV/i. 

The FCS is the probability distribution function P(Q) 
for the charges Q = {Qlt, Qrt, Qlb, Qrb) to 
be transferred through the respective junctions dur- 
ing a (long) waiting time r, thereby allowing to com- 
pute not only non-linear I-V, but also current noise 
and higher order cumulants. This information is en- 
coded in the CGF x(A) = Eq e**^^^(Q) where A = 
(Alt, Aljt, ^lb, ^rb) are the measuring fields. The 
cumulants (irreducible moments) are then found from 
the respective derivatives of Inx(A).— In order to ob- 
tain the CGF we adopt the approach of modifying the 
Hamiltonian by introducing a time-dependent counting 
field and relate xi^) to the Keldysh Green's functions of 
the systemJi Such procedure allows for the calculation 
of the FCS for arbitrarily given parameters of ([2|) and of 
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the tunneling amplitudes in (Iljji^ The determination of 
the currents and its cumulants in this hybrid structure, 
however, represents an essentially self-consistent prob- 
lem, where the external currents depend on the electro- 
chemical potentials of the two layers and on the excitonic 
order parameter, which in turn adjust to ensure charge 
conservation and no inter-layer tunneling, thereby affect- 
ing the external currents themselves. In order to pro- 
ceed, some assumptions are thus necessary. In view of 
possible implementations with graphene, we shall con- 
sider a linear Dirac cone spectrum Hq for the layers, 
oppositely shifted by e. g. two external gates ±Vg, so 
that Ht/b = -ffo T eVg - Hec,t/b, where Hec\t/b are 
the electrochemical potentials. It is sensible to focus 
on the incoherent tunneling regime, £,EcA<t> ^ ^ where 
Ifj, is the dephasing length. The condition £^ec *C I 
also implies that self-consistency effects on the space de- 
pendence of |A(a;)| are negligible 1^*^^ A space-dependent 
phase arg(A(a;)) ~ qx, on the other hand, although es- 
sential to ensure that the EC carries counterflowing cur- 
rents in the bulk of the bilayer, is not necessary for eval- 
uating the currents in the leads, which are of interest 
hereJ^ In contrast, self-consistency of the electrochem- 
ical potentials Hec,t/b of the two layers is crucial to 
ensure current conservation in each layer— 



{Ilt) — {Irt), 



{I lb) = {Irb). 



(3) 



Under these assumptions, we have obtained the com- 
plete analytical expression for the CGF for all parame- 
ter regimes. Such expression, which has been used for 
our numerical evaluation, is quite lengthy and we do 
not report it here. Nevertheless, all relevant ingredi- 
ents of the CGF already appear in the limits of small 
bias {fJ-L,a,lJ-R,a < A) and large bias {pL,a,^J-R,a > A), 
where the expression of the CGF greatly simplifies, and 
acquires the following form on the left leads 



Inxli =2r 



+ 



■' \a- = T, B 

ln{l + TA(aj) [O 



,g»A_f,^g->A, 



IN /I ^']\ n f ^-"^I^^H^tI, \i^B 

■ l)nLs(l -"-lt)! 1^ ' 



l)niT(l-niB) + (e'^^«e-'^ 



, (4) 



where the transmission coefficients arc given by T^{uj) ~ 
4fz,,/(l + TLaf and T^(w) = 4f^/(l + Ta)^- The 
effective transparencies are parametrised by the EC 
density of states as T^a = ^Lcr\^a\/ \/i^a ~ A^ and 



Ta = TLTTLBAyyA^ - wVA2 - w|], where Fl. = 
'^'^ PoLaPoElLa/'^- The functions /t and /s denote Fermi 
distributions for the quasiparticles in the separate lay- 
ers and cot.b = w — Pec,t/b- The first line of Eq. (|4]) 
describes the supra-gap contribution, which is only due 
to single electron transport and is characterized by the 
normal transmission coefficient T^. In contrast, the sec- 



I 

ond line describes the sub-gap contribution due to the 
phenomenon of excitonic Andreev reflection^, consisting 
of an electron and a hole (traveling in different layers), 
which enter or leave coherently the bilayer in order for 
an excitonic pair to be transfered along the bulk of the 
system. 

The expression for the currents in the left leads is 
{Ilct) = — ir^-'^ dhix/dXLa- Expressions for the rhs are 
obtained by replacing Alt — ^ — Aat, Alb -^ —^rb in 
Eq. (U]). Imposing the self-consistency condition ([3]) de- 
termines ^j-ecit, and one obtains the final results for the 
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FIG. 2. (color online) (a) current in the top (red) and bottom 
(blue) layer as a function of Vr, for a fixed value Vb = A/2. 
Flt = r_Ls = 0.171 corresponds to transmission of 0.5 for the 
uncoupled system, T = 0.01 A. For |Vt|, |Vb| < 2 A the bi- 
layer exhibits counterpropagating currents, exciton blockade 
occurrs at Vb = Vt- (b) differential conductance dlT/dVr 
(red) shows a resonance peak and tends to the typical value 
for a quantum point contact whereas the transconductance 
dJs/dVr shows a resonance peak before vanishing at larger 
bias. 



two currents !„ = {I La) = {Iro)- For simplicity, we 
consider the symmetric junction case, T lt = ^rt and 
Tlb = TflB, and symmetrically applied biases ^l^t = 
-M-RT = -Vt/2, i^iLB = -fJ-RB = -Vb/2. In this case 
Eq. ([3]) is always fulfilled for hect = I^ecb ~ 0. 
The average currents are plotted in Fig. [5] (a) as a func- 
tion of the top layer bias Vr, for a fixed value of the 
bottom layer bias Vb- As one can see, because of the 
EC coupling, both It and 1b change, even when varying 
Vt only. In particular, for |Vr|, |Vb| < 2A, one observes 
It = —Ib-, a signature that in the sub-gap regime trans- 
port can only occur via excitonic counterpropagating cur- 
rents in the bulk of the layers, which are transformed into 
electron and hole currents in the leads through excitonic 
Andreev reflections i^ Notice that for the value Vt = Vb 
a current locking occurs [Jt = Ib ^ 0), because the 
EC cannot sustain currents driven by equally applied bi- 
ases (exciton blockade). At Vr = 2 A excitonic pairs 
start to break up and the resulting electrons/holes get 
excited above the gap. This is clearly shown in Fig. [2] 
(b) , where the positive conductance exhibits a resonance 
peak, whereas the negative transconductance abruptly 
changes sign. At higher voltage values the EC plays a 
minor role, so that the conductance tends to the value of 
the case A = 0, and the transconductance vanishes, indi- 
cating that transport in the bottom layer is independent 
of the voltage applied to the top layer. 

The current correlators are defined as {{I Lai La')} = 
{laala'a') — {laa) {la' a') , and obtained from Eq. (m as 
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equilibrium (Vt — Vb — 0) we obtain the customary 



This result indicates that the two electrons involved in 
an excitonic Andreev reflection dwell in separate layers 
so that only the conductance of a single layer enters the 
Johnson-Nyquist noise. The second equality in ^ is 
the generalized Johnson-Nyquist relationship obtained in 
[2lj for a floating superconductor, which is here obtained 
without any use of Langevin forces. 
As we see from Fig. [3] also in nonequilibrium one al- 
ways observes a negative cross-correlation. This is differ- 
ent from the case of a superconductor contacted to two 
normal electrodes, where one observes a positive cross- 
correlation of the two currents in the normal leads via 
crossed Andreev rcflcctioni^ The reason for this crucial 
difference is the fact that in the EC case we probe the 
correlation of electrons and holes rather then correlations 
of electron pairs as in the case of superconductors. We 
call the two voltage bias situations Vt = ±Vb parallel 
and antiparallel configurations. In the first configuration, 
where the average current is vanishing, the noise and the 
cross-correlation are shown in Fig. [3Ja). In the sub-gap 
regime, up to thermal fluctuation effects, both the noise 
and the cross-correlation vanish. This is because the in- 
coming electrons and holes are always reflected back into 
the same lead they are injected from, since no exciton 
can penetrate inside the EC. Notice that this effect is 
essentially independent of the interface transmission, as 
shown in the three curves of Fig.[3][a). Indeed for the par- 
allel bias configuration in the subgap regime the EC gap 
effectively plays the role of a large barrier. In contrast, 
when \V\ > 2 A, quasiparticles can be excited above the 
gap, and the noise in each lead starts deviating from zero, 
eventually increasing linearly with V . 

In the antiparallel configuration, where the average 
currents in the layers flow in the opposite directions, the 
noise and the cross correlation are shown in Fig. ^h) 
as a function of V = Vt — —Vb- In the supra-gap 
regime the behavior is qualitatively similar to the parallel 
conflguration, so that the noise increases and the cross- 
correlations saturates to a value determined by the quasi- 
particle mixed character above the gap. However, differ- 
ences with respect to the parallel configuration emerge in 
the sub-gap regime, where both noise and cross correla- 
tion are now non-vanishing, and depend on the interface 
transmission. 

In realistic implementations of the proposed setup, a 
small single-particle inter-layer tunneling tab, here as- 
sumed vanishing, is expected to be present. The role 
of such a term has been discussed in the literaturo^ir— , 
and depends on the specific implementation. In quantum 
Hall bilayers with GaAs wells, the role of such term is to 
lift the degeneracy of the EC phase and to introduce a fi- 
nite critical current Ic,— which depends on the layer area 
and on |tabP- In typical drag setups, the total fiowing 
currents arc usually higher than Ic^ Inter-layer current 
is then negligible because the incoherent tunneling resis- 
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FIG. 3. (color online) (a) The parallel bias configuration (ex- 
citon blockade). Noise and cross-correlations are shown as a 
function of Vt = Vb = V, for the three values of transmis- 
sion for the uncoupled system of 0.3 (black), 0.5 (blue), 0.7 
(red) and T = O.OIA. (b) The antiparallel bias configuration 
(counterpropagating currents). Noise and cross correlations 
are shown as a function oi Vt ~ —Vb = V , again for the three 
values of contact transparency. 




FIG. 4. (color online) Transformation of voltages using the 
EC. We show Vb due to the current in the bottom layer 
coupled to a quantum point contact as a function of the 
applied voltage Vt for different transparencies of the QPC, 
Flt = 0.2 — Trt — Tlb = FflB and the QPC transparency 
is varied from Ti = 0.1 (blue curve), Ti = 0.3 (red) to Ti = 0.7 
(black). The inset shows the sketch of the experimental setup. 



tance is much higher than the in-plane resistance, so that 
the results presented here remain valid in the regime of 
currents Ic ^ |/|. We also checked that a long but fi- 
nite exciton lifetime leading to a small imaginary part of 
A does not alter our results significantly. In the case of 
graphene bilayers, the scenario can be richer, depending 



on the nature of the insulating barrier and on the rotation 
of the graphene layers with respect to each other. 

One possible application of the setup could be as a 
nanoscale voltage transformer. Ideal voltage transform- 
ers, such as inductors, require that the transformation 
coefficient depends weakly on the load characteristics, 
and that energy losses are minimal. However, imple- 
mentation of on-chip silicon based inductors turns out 
to be very difficulli^, so that the seek for trade-off solu- 
tions in nanotransformers represents a great challenge in 
modern electronics. The contact configuration proposed 
in Q suggests that the setup can be used as a voltage 
transformer at the nanoscale, contacting the EC bottom 
layer to another mesoscopic system with a known I-V 
characteristics, such as a quantum point contact (QPC), 
where Iqpc = GqTiVb, with Ti being the contact trans- 
parency (see inset of Fig. 2]). At low temperatures and 
for \Vt, b\ < 2A, only excitonic Andreev reflections con- 
tribute to transport and determine (Ilt) as a function 
of Vt- On the other hand (Ilt) = —{Ilb), and current 
conservation implies that —{I lb) = Iqpc- This leads to 
the transformation 



{Ilt){Vt)^GoTiVb- 



(6) 



A typical example for the voltage interrelation is shown 
in Fig. m In this case the transformation is controlled 
by Ti, which is typically tunable. Other realizations may 
involve transistors based on coupling to internal degrees 
of freedom. Importantly, nanotransformers based on the 
dissipationless EC counterflowing currents may help min- 
imizing heat and noise production. 

To conclude, we have calculated the FCS for an 
EC contacted to four metallic leads. We have shown 
how counterpropagating currents and the generalized 
Johnson-Nyquist relation directly follow from the cumu- 
lant generating function. Using this approach we ana- 
lyzed noise driven by excitonic Andreev reflections. Al- 
though the effective model is quadratic in fermion fields it 
correctly describes the non-trivial multi-particle exciton 
bound states. The lowest cumulants of charge transport 
resemble the free electron result, nonetheless fully ac- 
counting for the highly non-trivial drag effects as well 
as energy-dependence of the effective transmission co- 
efficients. We also showed how excitons can be used 
for transforming current on the nanoscale. The authors 
would like to thank S. Maier, R. Fazio, A.H. MacDonald, 
M. Polini, and F. Taddei for many interesting discus- 
sions. The financial support was provided by the DFG 
under grant No. KO-2235/3, and 'Enable fund' of the 
University of Heidelberg. 
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